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Abstract 

High energy neutrinos play a very important role for the understanding of the 
origin and propagation of ultra high energy cosmic rays (UHECR). They can be 
produced as a consequence of the hadronic interactions suffered by the cosmic rays 
in the acceleration regions, as by products of the propagation of the UHECR in 
the radiation background and as a main product of the decay of super heavy relic 
particles. A new era of very large exposure space observatories, of which the JEM- 
EUSO mission is a prime example, is on the horizon which opens the possibility of 
neutrino detection in the highest energy region of the spectrum. In the present work 
we use a combination of the PYTHIA interaction code with the CONEX shower 
simulation package in order to produce fast one-dimensional simulations of neutrino 
initiated showers in air. We make a detail study of the structure of the corresponding 
longitudinal profiles, but focus our physical analysis mainly on the development of 
showers at mid and high altitudes, where they can be an interesting target for space 
fluorescence observatories. 
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1 Introduction 



The high [Ei^ < 10^^ eV) and ultra high [E^, > 10^^ eV) energy neutrino 
fluxes carry very important astrophysical information. In particular, neutrinos 
arriving to the Earth can be originated in very distant sources because they 
travel through the universe without interacting. A high energy neutrino flux 
is expected as a by-product of the interactions of cosmic ray hadrons at the 
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sources (see e.g. [1]). They can also be produced during the propagation of 
cosmic rays in the intcrgalactic medium [2,3] and as the main product of the 
decay of superheavy rehc particles [8,9]. 

There are essentially two different strategies to detect neutrinos in cosmic 
rays detectors. The first consist in observing the development of horizontal air 
showers produced by the interactions of electron neutrinos with nucleons of the 
molecules of the Earth atmosphere [7] , and the second one consist in observing 
the showers produced by the decay of taus generated by the interaction of tau 
neutrinos propagating through the interior of the Earth [4,5,6]. 

Space observatories play a very important role in neutrino detection, in par- 
ticular, JEM-EUSO [10] with its lO-*^^ tn of atmospheric target volume has the 
real possibility of observing ultra high energy neutrinos and make important 
contributions to the understanding of UHECR production and propagation 
[11,12]. Source distributions rapidly evolving with redshift would be particu- 
larly favorable by increasing the cosmogenic neutrino flux at highest energies 
[2]. A thorough understanding of neutrino deep inelastic scattering, as well 
as the evolution of longitudinal profiles of atmospheric neutrino showers, are 
extremely important in order to take advantage of the full potential of the 
experiment. Conversely, besides the obvious astrophysical value, the proper- 
ties of just a few observed showers can also give valuable information on the 
physics governing high energy neutrino-nucleon interactions. 

In this work we first study the neutrino-nucleon interactions at the highest 
energies for two different sets of parton distribution functions. We study in 
detail, the characteristics of horizontal air showers originated by the interac- 
tions of electron neutrinos in the Earth atmosphere. In particular, we consider 
horizontal showers initiated at different altitudes and very deep in the atmo- 
sphere, which, depending on the flux, will be detected by the upcoming orbital 
detectors like JEM-EUSO. We also study the detectability of such showers as 
a function of the altitude, for an ideal orbital detector similar to JEM-EUSO. 



2 Neutrino nucleon interaction 

High energy neutrinos that propagate in the Earth atmosphere can interact 
with protons and neutrons of the air molecules. There are two possible channels 
for this interaction, charged and neutral current, vi-\-N ^l-\-X and vi-\-N ^ 
1^1 + X, respectively. Here N is a nucleon (proton or neutron), i/i is a neutrino 
of family /, / is the corresponding lepton and X the hadronic part of the 
processes. At the level of the quark-parton model, the entire hadronic state 
of a deep inelastic scattering may be viewed as the fragmented product of 
a scattered quark and the proton remnant. The major uncertainty on the 
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differential cross section at the energies considered comes from the unknown 
behavior of the parton distribution functions (PDFs) at very small values of 
the parton momentum fraction x [13]. 

In this work, the simulation of the neutrino nucleon interaction is performed 
by using the PYTHIA code [14]. PYTHIA is an event generator, intended 
for high-energy processes with particular emphasis on the detailed simulation 
of quantum chromodynamics (QCD) parton showers and the fragmentation 
process. 

The parton shower approach was developed to take into account higher than 
first order QCD effects. It has the advantage that arbitrarily high orders in 
the strong coupling constant can be simulated, but only in the leading or- 
der approximation, as opposed to the exact treatment in fixed order matrix 
element. Higher order effects are important at high energies where multiple 
parton emission can give rise to multijet events as well as affect the internal 
properties, such as hardness and width, of a jet (see Ref. [15]). 

QCD perturbation theory, formulated in terms of quarks and gluons, is valid at 
short distances. At long distances, it becomes strongly interacting and pertur- 
bation theory breaks down. As mentioned before, in this confinement regime, 
the colored partons are transformed into colorless hadrons. The fragmentation 
process has yet to be understood from first principles, starting from the QCD 
Lagrangian. As a consequence, a number of different phenomenological mod- 
els have been developed to describe this effect. In PYTHIA the fragmentation 
process is done by using the so-called Lund string model [16]. 

A typical high energy event has the following structure [14]: 

(1) At the beginning of the simulation two incident beams are coming in 
towards each other. Each particle is characterized by a set of PDFs which 
determines the fraction of momentum taken by each parton. 

(2) A collision between two partons, one from each beam, gives the hard 
process of interest. A collision implies accelerated color (and often elec- 
tromagnetic charges), therefore, bremstrahlung can occur. The colliding 
partons start off a sequence of branchings (such q ^ qg, g ^ gg, g — >■ qq, 
etcetera) which build up an initial-state shower. 

(3) Also the outgoing partons may branch to build up a final-state shower. 

(4) At this stage just one parton of each incident beam is taken out to undergo 
a hard collision. The beam particles are made up of a multitude of further 
partons, then, a beam remnant is left behind. This remnant may have an 
internal structure and a net color charge. 

(5) The QCD confinement mechanism ensure that the outgoing quarks and 
gluons are not observable, but instead they fragment to color neutral 
hadrons. 
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Fig. 1. Average energy fraction taken by the electron in a proton-electron neutrino 
charge current interaction obtained from PYTHIA with CTEQ6 (circles) and GJR08 
(squares) sets of PDFs. 

(6) Many of those primary hadrons are unstable and decay further at various 
time scales. 



The default configuration of PYTHIA is used in the present simulations. 
Also, the parton distribution library LHAPDF [17] is linked with PYTHIA 
to use different extrapolations of the PDFs. In order to study the influence of 
the PDFs on the electron neutrino showers, two different sets are considered: 
CTEQ6 [18], the most commonly used in the literature (at the highest ener- 
gies) and GJR08 [19]. Fig. 1 shows the energy fraction carried by the electron 
as a function of the incident neutrino energy for the charge current interac- 
tion of an electron neutrino with a proton for both sets of PDFs considered. In 
both cases the energy fraction increases steadily with the incoming neutrino 
energy. The difference between both PDFs increases up to a maximum of a 
few percent at the highest energies. Fig. 2 shows the distribution of the energy 
fraction taken by the electron for a neutrino of = 10^° eV interacting with 
a proton corresponding to CTEQ6. On average, approximately 82% of the 
neutrino energy is taken by the electron. 

Besides the leading particle, different types of secondaries are generated as a 
result of the interaction. In particular, in this work we are interested in the 
ones recognized by CONEX [20] code which is used to simulate the neutrino 
showers (see section 3). The few final state particles that are not treated yet 
by CONEX are mapped into their decay products [21]. Fig. 3 shows the energy 
fraction taken by the most relevant particles recognized by CONEX for three 
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Fig. 2. Distribution of the energy fraction taken by the electron in a electron neutrino 
proton interaction. The neutrino energy is E^, = lO'^^ eV and PYTHIA with CTEQ6 
sets of PDFs are used for the simulations. The number of simulated events is 1000. 
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Fig. 3. Average energy fraction for the most important particles recognized by 
CONEX produced in a charge current interaction of an electron neutrino with a 
proton for three different neutrino energies. 

different electron neutrino energies, obtained as a result of the charge current 
interaction with a proton. It can be seen that the smaller the neutrino energy 
the larger the energy fraction taken by the secondary particles. 

Note that the results obtained for electron neutrinos are also valid for electron 
antineutrinos because at the energies considered in this work the cross-sections 
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are approximately the same (see for instance Ref. [22]). 



3 Neutrino showers 

There are several method to study the development of the extensive air show- 
ers. The most common one is the Monte Carlo (MC) simulation of the interac- 
tions suffered by the primary and the secondary particles of the showers with 
the air molecules. The MC programs most used in the literature are COR- 
SIKA [23] and AIRES [24]. At the highest energies, the number of particles 
of the cascades is extremely large, e.g. ~ 10^^ particles at the maximum of a 
10^° cV proton shower. The computing time required to simulate the interac- 
tions of all particles is so large that sampling algorithms, so-called "thinning" , 
are used [25,26]. They are such that just a small fraction of the particles are 
propagated and a weight factor is associated to each of them. In these MC 
programs, both, the lateral and the longitudinal development of the showers 
are simulated. 

The hybrid approach is a possible alternative to speed up the simulation of 
the showers without losing accuracy. It consists in the detailed MC simula- 
tion of the cascade for particles with energy above a given threshold and a 
description of the low energy sub-showers based on the numerical solution of 
the corresponding cascade equations. 

As mentioned before the CONEX program is used to generate the electron 
neutrino showers. CONEX is a one-dimensional hybrid program which can 
be used to simulate the longitudinal profile of the showers very fast and very 
accurately. The accuracy of the calculation is supported by comparisons done 
between CONEX simulations and CORSIKA ones (see Ref. [20] for details). 
Note that the lateral development of the showers, not included in CONEX, is 
not important for the observation from the space, the showers are very well 
approximated by a point for the distances involved in this kind of techniques. 

The major uncertainty for the shower simulations comes from the unknown 
of the hadronic interactions at the highest energies. There are several models 
that extrapolate the accelerator data to the energies of the cosmic rays, one of 
those is QGSJET-11 [27] which is the one used in this work. The development 
of the shower strongly depends on the assumed hadronic interaction model. 
In particular, the most important observables of the showers, hke the position 
of the maximum and the muon content, are quite sensitive to the hadronic 
interaction model considered (see Ref. [28,29] for details). 

The particles produced in a neutrino-nucleon interaction are injected in CONEX 
with QGSJET-II producing extensive air showers. The energy thresholds used 
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Fig. 4. Horizontal electron neutrino showers at four different altitudes of Ei, = 10 
eV, injected at a point contained on the vertical axis of a nadir-pointing orbital 
detector. = 15 profiles are shown for every value of altitude considered. 

in the present simulations are the ones suggested by the authors of the pro- 
gram, i.e. default configuration. 

Because the mean free path of neutrinos propagating in the atmosphere is 
very large, they can interact very deeply, after traversing a large amount of 
matter. An orbital detector like JEM-EUSO can also detect horizontal showers 
that do not hit the ground. In particular, horizontal neutrinos can interact 
at higher altitudes producing a shower observable by the detector. Fig. 4 
shows the energy deposit as a function of X — Xq, where Xq corresponds to 
the atmospheric depth of the injection point, for horizontal electron neutrino 
showers of E^, = 10^° eV. The injection points are such that the trajectory 
of the showers starts at the vertical axis of a nadir-pointing orbital telescope 
(like JEM-EUSO in nadir mode, see section 4) and at different altitudes. Note 
that, because of the very large mean free path of neutrinos in the atmosphere, 
a subset of all possible neutrino interactions inside the field of view (FOV) of 
the detector can be modeled by this kind of trajectories. Also note that for an 
horizontal proton shower oi E = 10^" eV the maximum is reached at X = 890 
g cm~^, with an energy deposit of ~ 1.65 x 10® GeV g^^ cm^. 

As seen from fig. 4, the profiles corresponding to smaller altitudes are very 
broad profiles which may present several peaks and large fluctuations. This 
behavior is due to the Landau Pomeranchuk Migdal (LPM) effect, which is 
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Fig. 5. Mean value and one sigma regions of the longitudinal profile corresponding 
to horizontal electron neutrino showers of Ey = IQ^^ eV for different altitudes. The 
interaction point is contained on the vertical axis of the JEM-EUSO telescope. 

very important inside dense regions of the atmosphere and for electromagnetic 
particles, electrons in this case, which take about 80% of the parent neutrino 
energy. Fig. 5 shows the mean value and one sigma regions of the longitudinal 
profiles for the same geometry and injection point considered before. It can 
be seen that as the altitude increases the fluctuations are reduced and, on 
average, the profiles become thinner. This is due to the fact that the LPM 
effect become progressively less important with increasing altitude because of 
the decrease in atmospheric density. Note that just showers up to 20 km of 
altitude are considered because, at higher altitudes, the grammage of the FOV 
like the one of JEM-EUSO, 60°, is not enough to contain the whole profiles. 

As already mentioned, at smaller altitude the longitudinal profiles present a 
complicated structure (see Fig. 4). In particular, the showers present multiple 
peaks. Fig. 6 shows the distributions functions of the position of each maxi- 
mum. X^^^ corresponds to the position of the first maximum counted from 
the start of the shower, X^^^. is the second one and so on. The maxima are 
obtained by searching for the points of the longitudinal profile which have six 
consecutive neighbors, three with atmospheric depth smaller and three with 
atmospheric depth larger. These six points are such that the deposited energy 
is smaller than the one of the candidate. When a candidate is found, the coor- 
dinates of the maximum are obtained by using a parabola that fit the central 
point and the two nearest neighbors, one at the right and the other at the left. 
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Fig. 6. Distribution of the position of the maxima for electron neutrino showers of 
Ey = 10^0 eV at sea level. 
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Fig. 7. X^^^ — Xq as a function of the fractional energy taken by the electron, for 
neutrino showers of Ei, = 10^^ eV at sea level. 



Note that the distribution function of X^^^—Xq is bi- valued and its first peak 
is located at ~ 800 g cm~^, while the second one is centered at ~ 1500 g cm~^. 
The first peak corresponds to the development of the hadronic component 
of the electron neutrino cascade, whereas the second one mainly reflects the 
electromagnetic portion of the shower. Fig. 7 shows X^^^ — Xq as a function of 
the energy fraction taken by the electron, it can be seen that, as the fractional 
energy of the electron decreases X^^^ — Xq of the events approaches to the 
values corresponding to the first peak. 

The interpretation of the X^^^ — Xq distribution can be confirmed by simu- 
lating the shower development removing the hadronic component produced in 
the neutrino nucleon interaction, i.e. considering just the generated electron. 
The distribution of X^^^ — Xq for the electron showers obtained applying the 
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same analysis as before do not present the first peak, as expected it just have 
one peak at ~ 1500 g cm~^ (see Fig. A.l of appendix A). 

The simulations show that ~ 65% of the showers that have a second maximum 
have X^^^ — Xq < 1000 g cm^^. Therefore, as mentioned before, the distri- 
bution of X^^^ — Xq has to do with the electromagnetic part of the showers. 
When the first maximum is in the electromagnetic region of the distribution, 
the second one is originated by the fluctuations due to the LPM effect. This 
is supported by the X^^^ — Xq distribution of the electron showers discussed 
in the appendix A, it is originated by the LPM fluctuations and, as expected, 
corresponds to larger values of X^^^^ — Xq (see Fig. A.l). The third maximum 
i-^max) is always originated by the fluctuations caused by the LPM effect. 

This complicated structure simplifies as the altitude increases. Fig. 8 shows 
the probability of finding a profile with a given number of peaks, A^x»„„^- As 
the altitude increases the probability of finding a shower with more than one 
peak goes to zero. 

The first portion of the cascades is, in general, dominated by the hadronic 
component. Therefore, the way in which the primary neutrino energy is dis- 
tributed can be assessed by. 
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Fig. 9. Distribution function of parameter Fen (see Eq. (1)) for CTEQ6, GJR08 
models and also for an artificial modification of the CTEQ6 prediction such that 
the leading particle takes 70% of the neutrino energy. Horizontal showers of 10^*^ 
eV injected at sea level at the axis of the FOV are considered. 

where dE/dX is the energy deposition, X^. = 1000 g cm~^ is a characteristic 
depth that roughly separates the hadronic-dominated from the electromagnetic- 
dominated portions of the shower (see Fig. 6), and Xum is the maximum at- 
mospheric depth visible by a given orbital detector. Fen is calculated only for 
those showers that present a first maximum at a depth smaller than Xe- Fig. 
9 shows the Fen distributions obtained for CTEQ6, GJR08 and for a modifi- 
cation of the CTEQ6 results in which the leading particle takes 70% of the 
neutrino energy, corresponding to horizontal lO^'' eV electron neutrinos in- 
jected at sea level at the axis of the FOV. It can be seen that Fen is correlated 
with the energy taken by the leading particle. The smallest value of the av- 
erage of Fen corresponds to CTEQ6 because the leading particle takes about 
82% of the neutrino energy while, in the case of GJR08, it takes ~ 79%. The 
differences between the results obtained for these models are very small. In the 
case in which the energy taken by the leading particle is artificially reduced 
to 70% (hashed histogram), the mean value of Fen increases. Note that, in the 
latter case, the energy extracted from the leading particle in order to reduce 
its average energy, is redistributed among the other daughter particles in such 
a way that the ratio between the energy taken by a given particle and the total 
energy taken by all secondaries, excluding the leading particle, is constant. 

Figures 7 and 9 show the influence of the energy taken by the hadronic com- 
ponent of the neutrino interaction on the hadronic character of the beginning 
of electron neutrino showers. In fact, the position of X^^^ — Xq as well as the 
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energy deposition at the beginning part of the showers are correlated with the 
amount of energy that goes to this component. 

The sensitivity of Fen to the hadron component of the showers depends on the 
altitude. Showers injected at sea level in regions of high density are dominated 
by the LPM effect allowing a clearer separation between the hadronic and 
electromagnetic portions of the cascades. 

The parameter Fe„ can be very useful to understand neutrino interactions 
with atmospheric nuclei and, in particular, to estimate the energy fraction 
taken by the leading particle. Any practical application, however, will depend 
on the actual event rate. 



4 Detectability for an ideal orbital detector 

As mentioned before, a new generation of orbital detectors are planned. Such 
instruments are designed to detect the fluorescence light emitted by the inter- 
action of the charge particles of the air showers initiated by the primary cosmic 
rays. These new instruments, of which JEM-EUSO is the pioneer, watch the 
atmosphere from the space, having a huge effective area. In particular, JEM- 
EUSO will watch the atmosphere in two different ways, nadir and tilted modes. 
The telescope is formed by three double-sided Fresnel lenses made of plastic 
material and an aspherically curved focal surface composed by multi-anode 
photomultipliers to detect the incident photons [30]. 

A simplified simulation of an orbital detector with characteristics similar to 
JEM-EUSO is developed, in order to study the detectability of the electron 
neutrino showers. The simulation of the fluorescence hght emitted by the 
charged particles of the cascades is done by using the photon yield of Ref. 
[31]. The fluorescence light is propagated from the shower axis to the tele- 
scope considering the Rayleigh and Mie attenuation [32]. An ideal optical 
system is considered, an optical transmission of 100% is assumed independent 
of the angle and wavelength of the incident hght. Regarding the PMTs a com- 
bined collection efficiency and quantum efficiency of ep^j = 0.27 is considered 
[33]. The Cherenkov photons are not included in the present simulations, for 
the case of horizontal showers this component is less important because such 
photons are very coUimated with the shower axis and just the scattered ones 
can contribute to the hght collected by the telescope. Also the attenuation 
of the fluorescence light due to the ozone absorption is not included in the 
simulation. 

Table 4 shows the parameters used in the simulation of the ideal orbital de- 
tector (see Ref. [34]). 
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Parameter 



Value 



Height of the orbit 



400 km 



Aperture diameter {Da) 



2.5 m 



Field of View 



±30° 



Wavelength range 
Gate Time Unit (GTU) 
Number of pixels (Npix) 



330 - 400 /xm 



2500 ns 



2 X 10^ 



Pixel size 



0.1 



o 



Table 1 



Parameters used in the simulation of an ideal orbital detector similar to JEM-EUSO. 

The main source of background light is the atmospheric nightglow, the mean 
value assumed in this work is Bq = 500 photons m~^ ns~^ sr~^ [30]. The rate 
of background photons per pixel is given by, 



where A — nD^/A and 7/o„ = 30° (7 is used to note the angle between the 
incident hght and the axis of the telescope) . 

Figure 10 shows the number of photons that produce signal in the PMTs 
as a function of the arrival time, for two different events corresponding to 
horizontal electron neutrino showers injected at sea level in the center of the 

FOV of the telescope. The temporal width of each bin is equal to a GTU. 
The black lines, N{k), correspond to 3 (solid line) and 5 (dashed line) sigmas 
above the background in the corresponding time bins. Assuming Poissonian 



statistics and for k sigmas above the mean in the ith bin, Ni{k) = k J Nf, with 



N^^RhX GTU X A7i/0.1°, where A7i/0.1° = {-/{U + AU) --/{ti)) /0.1° is the 



fraction of the pixel corresponding to the ith time bin. The separation between 
two consecutive dotted vertical lines indicates the time interval corresponding 
to an angular separation of 0.5°, i.e. five pixels. 

The event on the left panel of Fig. 10 presents just one peak, it has several 
pixels with the number of photons well above the background, 12 above 3 a 
and 5 a. The event on the right panel of the same figure has two peaks, in this 
case the photons are distributed in a wider interval of time due to the shape 
of the longitudinal profile. Also in this case there are several pixels with the 
number of photons above the background, 23 above 3 a and 18 above 5 a. 

As the altitude increases the distribution of photons changes. Fig. 11 shows 
two events corresponding to /iq = 5 km, it can be seen that the number of 
photons that produce signal in the PMTs increases because the distance of the 
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Fig. 10. Time distribution of photons that produce signal in the PMTs for two 
electron neutrino horizontal showers of Ei, = 10^^ eV, injected at sea level in the 
center of the FOV of an orbital telescope. The set of PDF used is CTEQ6. 
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Fig. 11. Time distribution of photons that produce signal in the PMTs for two 
electron neutrino horizontal showers of Ei, = W"^^ eV, injected at an altitude of 
ho = 5 km in the center of the FOV of an orbital telescope. The set of PDF used is 
CTEQ6. 



shower to the telescope is smaller, increasing the solid angle of the photons that 
reach the telescope and decreasing the attenuation in the atmosphere. Also 
the angular width of the distribution increases due to geometrical effects. As 
a result, the number of pixels with the number of photons larger than a given 
level of background increases. As expected, the detectability of the horizontal 
showers improves for larger altitudes. For the event on the left panel of the 
figure there are 22 pixels above 3 a and 20 pixels above 5 a and for the event 
on the right panel there are 38 pixels above 3 a and 33 pixels above 5 a. 

Figure 12 shows the distribution of photons for showers injected at /io = 10 
km and Hq = 20 km. Again, the number of pixels with signal well above a 
given background level increases with altitude in such a way that, for ho = 10 
km (left panel) there are 44 pixels above 3 a and 39 pixels above 5 a. For the 
event corresponding to Hq = 20 km (right panel) there are 187 pixels above 
3 a and 171 pixels above 5 a. 
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Fig. 12. Time distribution of photons that produce signal in the PMTs for two 
electron neutrino horizontal showers of Ei, = 10^^ eV, injected at ho = 10 km (left 
panel) and /iq = 20 km (right panel) in the center of the FOV of an orbital telescope. 
The set of PDF used is CTEQ6. 

5 Proton and neutrino showers 



The interaction length for protons is Apr(lO^'^eV) ~ 36 g cm^^ (for Sibyll 2.1 
[36]) and for neutrinos is A,^(10^'^eV) ~ 3.2 x 10'' g cm~^. The survival prob- 
ability of an horizontal proton that reaches the Earth surface at the vertical 
axis of the FOV is ~ exp(— 1000), whereas the corresponding probability for 
a neutrino is ~ exp(— 0.001). Therefore, despite the fact that horizontal neu- 
trino and proton showers have very different observational characteristics, it 
is very unlikely to observe a proton interacting so deep in the atmosphere. 

Nevertheless, for a given proton and a neutrino fluxes, there exists a particular 
slant depth for which the proton and neutrino events have the same rate. For 
any particle type, the probability of interacting in the interval [X, X + AX] 
is given by, 

P^E; X, AX) = exp{-X/X{E)) [1 - exp(-AX/A(E))] , (3) 

where X{E) is the interaction length at a given energy. Therefore, solving the 
equation 0p,(E) PZiE] Xq, AX) = (f)^{E) P^^t{E- Xq, AX), where 0p,(E) and 
4>iy{E) are the proton and neutrino fluxes, the slant depth at which protons 
and neutrinos can be detected with the same rate is, 



j^^^^ _ Kr{E)\u{E) 



Xpr{E) — XuiE) 



( ME) \ / l-exp(-AX/A,(E)) Y 
[M{E))^ [l~exp{-AX/X,r{E))) 



(4) 



The function X{E) is obtained by using: AX = Xpr{E), the proton-air cross 
section of Sibyll 2.1, the charged current neutrino cross section a%'^{E) = 
6.04 X 10-3^ (E/GeV)°-3^^ cm^ [37], the Waxman-Bachall upper limit for the 
neutrino flux [38] and a power law flt of the Auger spectrum [39]: (ppr{E) = 
3.64 X 10^* (E/eV)-3-733 ni-^ sr^^ eV~^ Therefore, for E = 10^° eV the 
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Fig. 13. Proton and electron neutrino air showers of 10 eV, 9 = 60° injected at 
360 g cm~^. 

event rate for protons and neutrinos are of the same order of magnitude for 
X = 360 g cm-2 (where ^^.(lO^OeV) = 2 x 10"^^ m'^ s"^ sr"^ eV"^ and 
0pr(lO^°eV) = 8 X 10-3^ m-2 s"^ sr'^ eV"^). The latter means that, under 
the assumptions of the present calculation, protons can act as a background 
for neutrino identification in inclined events. Note that for values low enough 
of the neutrino flux, which could be the case, there is no atmospheric depth 
for which both rates are the same. 

Fig. 13 shows the average profile and the one sigma region for proton and 
electron neutrino induced air showers of zenith angle 6 = 60° and primary 
energy lO^'' eV injected at 360 g cm~^. Although, the event rates of this kind 
of showers are similar by assumption, for both protons and neutrinos, the 
profiles, however, are quite different making possible a good discrimination. 
More specifically, it can be seen that the neutrino showers develop deeper in 
the atmosphere and present larger fluctuations, due to the LPM effect. 



6 Conclusions 

Neutrino detection is of great importance for the understanding of several 
astrophysical process and, in particular, the origin and propagation of the 
highest energy cosmic rays. Orbital detectors like JEM-EUSO are technically 
capable of observing neutrino initiated cascades and to discriminate these 
cascades from hadronic ones. 
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Horizontal electron neutrino showers interacting very deep in the atmosphere 
are dominated by the LPM effect. The characteristics of such showers depend 
strongly on the density of the atmosphere and consequently on the altitude 
at which they develop. We show that horizontal electron neutrino showers 
interacting at sea level can present several peaks and large fluctuations. As 
the altitude increases the probability of finding a profile with more than one 
peak decreases. Also the fluctuations decrease with altitude. We find that the 
detectability of this kind of showers, in the context of an orbital telescope, 
improves with altitude and the multiple peak structure can be observed by 
this type of detectors. We also show that protons are a possible background for 
the neutrino identification depending on the unknown neutrino flux. However, 
even for optimistic values of such flux the longitudinal profiles are different 
enough to allow a good discrimination. 
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A Electron showers 

Electron showers are generated injecting in CONEX just the electron gener- 
ated in the electron neutrino nucleon interaction simulated with PYTHIA. 
The distributions of X^^^ — Xq are obtained by using the same method de- 
scribed in section 3. Figure A.l shows that, in this case, the distribution of 
the first maximum just has one peak at ~ 1500 g cm~^. 
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